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1.  INTRODUCTION 


The  main  aim  of  this  study  is  to  design  and  demonstrate  a  planar 
optical  waveguide  version  of  an  optical  isolator  employing  the  Faraday 
effect.  This  effect  produces  a  rotation  of  the  plane  of  polarization 
of  linearly  polarized  light  as  it  passes  through  a  bulk  magneto-optic 
medium.  The  rotation  is  right-handed  (say)  for  light  travelling  with 
the  magnetic  field  and  left-handed  (say)  for  light  against  the  field. 
This  non-reciprocal  nature  has  been  exploited  at  microwave  ( waveguide 
modes)  and  optical  (bulk-wave  modes)  frequencies  to  make  one-way 
transmission  devices  or  isolators^1'.  The  situation  is  complicated 
in  planar  optical  waveguides  because  circularly  polarized  propagation 
modes  (vhich  are  used  to  explain  the  optical  rotation  of  linearly 
polarized  light)  do  oot  normally  propagate  in  planar  waveguides. 

In  this  annual  technical  report  we  first  of  all  summarize  the 
essential  features  of  a  bulk-wave  optical  isolator,  indicating 
what  refinements  over  the  basic  configuration  are  required  before  a 
dielectric  waveguide  analog  can  be  conceived.  We  then  outline 
our  path  towards  a  practical  design  of  dielectric  waveguide  isolator. 
This  particular  design  has  been  given  a  thorough  theoretical  analysis, 
as  outlined  in  section  3,  and  a  practical  version  based  on  magnetic 
garnet  films  has  been  the  subject  of  our  experimental  program.  We 
give  design  data  in  section  4  and  at  the  end  of  the  report  give 
a  summary  of  our  achievements  to  date.  Certain  experimental  aspects 
of  the  program,  such  as  the  technique  for  determining  waveguide 
thickness  and  refractive  inF  ;  parameters,  should  find  application 
outside  the  scope  of  this  immediate  study,  and  these  techniques  are 
described  in  section  5  of  the  report. 

2.  BASIC  DESIGN  CONSIDERATIONS 


2.1.  Bulk-wave  Optical  Isolators 

Figure  1  illustrates  the  longitudinal  magneto-optic,  or 
Faraday  effect  which  is  exploited  in  the  design  of  optical  isolators 
and  circulators.  An  optical  beam  is  sent  through  a  magneto-optic 
crystal  parallel  to  an  applied  magnetic  field.  Linearly  polarized 
light  travelling  with  the  field  will  find  its  polarization  direction 
rotated  as  a  right  hand  screw  whereas  light  travelling  against 
the  field  is  rotated  in  the  opposite  sense  (a  left  hand  screw  in  our 
example).  This  rotation  is  similar  to  that  experienced  in  optically 
active  crystals  with  one  important  difference;  the  magnetic  medium 
experiences  right-  and  left-hand  rotation  whereas  the  handedness 
of  rotation  in  an  optically  active  medium  is  not  changed  by  reversing 
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the  propagation  direcLion.  If  the  path  length  in  the  magneto-optic 
crystal  is  such  that  a  45°  rotation  is  achieved  we  are  able  by  the 
simple  addition  of  two  linear  polarizers  at  suitable  azimuth 
orientations,  to  devise  "one-way"  transmission  characteristics  which 
constitute  an  optical  isolator.  Figure  2  shows  that  the  polarizers 
are  set  at  45°  to  each  other.  Vertically  polarized  light  travelling 
with  the  applied  magnetic  field  passes  through  the  input  polarizer, 
is  then  rotated  with  a  right-hand  screw  motion  to  45°  azimuth  where 
it  is  correctly  oriented  to  pass  through  the  second  polarizer.  If 
this  45°  polarized  light  is  sent  back  against  the  field  direction  it 
would  be  rotated  with  a  left-hand  screw  motion  to  a  90°  azimuth  to  be 
blocked  by  the  polarizer  at  0°  azimuth.  In  this  way  a  combination 
of  two  polarizers  set  at  45°  to  each  other  and  a  45°  Faraday  rotator 
will  act  as  an  optical  isolator,  permitting  transmission  in  one 
direction  only. 

Waveguide  modes  in  planar  dielectric  waveguides  have  their 
electric  vectors  vibrating  in  orthogonal  planes.  The  simple  bulk 
isolator  described  above  will  not  therefore  carry  over  to  a  waveguide 
version  because  the  input  and  output  polarizations  are  at  45°,  not 
0  or  90°.  A  second  active  element  must  be  included  so  that  the  output 
polarization  is  brought  back  to  0°  azimuth  as  illustrated  in  Fig.  3. 

A  45°  optically- active  rotator  with  a  left-hand  sense  of  rotation  is 
inserted  between  the  magnetic  rotator  and  the  second  polarizer. 

This  brings  light  vibrating  at  45°  azimuth  to  the  vertical  for  both 
directions  of  propagation.  Thus  v. '  can  see  intuitively  that  we  need 
two  active  media  to  make  an  optical  waveguide  isolator.  It  is  also 
implicit  in  those  bulk  wave  isolators  that  the  vertical  and  horizontal 
polarizations  travel  in-phase  and  the  observed  rotation  of  the  direc¬ 
tion  of  vibration  is  a  consequence  of  the  coupling  between  these 
polarizations.  This  phase -matching  requirement  carries  over  to  the 
optical  waveguide  cases  where  the  TE  and  TM  modes  must  travel  ir.  phase. 

2 . 2  Planar  Dielectric  Waveguide  Isolators 

(2) 

A  waveguide  version  of  Fig.  3  has  been  proposed  by  Wang  et  al 
in  which  magneto-optic  and  optically  active  (or  anisotropic)  crystals 
are  placed,  in  tandem,  on  top  of  a  dielectric  waveguide  as  shown  in 
Fig.  4.  The  waveguide  mode  fields  penetrate  into  these  regions  but 
decay  exponentially  n't  they  do  so.  It  is  possible  to  replace  the 
optically  active  crystal  with  an  anisotropic  one  in  optical  waveguide 
cases  because  TM  modes  have  a  component  of  electric  field  along  the 
propagation  direction  and  an  anisotropic  crystal  can  be  oriented  to 
provide  a  coupling  term  between  this  field  component  and  the  electric 
field  of  a  TE  mode. 

The  magnetic-optic  effect  only  provides  a  weak  coupling  between 
TE  and  TM  fields  and  since  in  Fig.  4  there  is  only  a  small  penetration 
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of  the  mode  field  into  the  top  layers  the  amount  of  mode  conversion 
attainable  would  be  quite  small.  This  in  turn  calls  for  a  requirement 
of  accurate  phase-matching  between  the  TE  and  TM  modes  over  long 
distances  (on  the  order  of  10  to  20  mms).  Mode  degeneracy  cannot  be 
achieved  in  a  dielectric  waveguide  mode  with  optically  isotropic 
materials  (see  Section  3.2)  and  so  the  waveguide  structure  supporting 
the  magneto-optic  and  anisotropic  crystals  would  itself  have  to  be 
made  from  anisotropic  materials. 

A  much  better  basic  design  is  achieved  if  the  essential  materials 
(i.e.  the  magneto-optic  and  anisotropic  layers)  are  regrouped  as 
shown  in  Fig.  5.  In  this  configuration  the  magneto-optic  film 
forms  the  high  index  guiding  layer  and  consequently  is  situated 
where  the  waveguide  modes  are  most  strongly  located.  Also  the 
anisotropic  layer,  which  is  placed  as  a  supeistrate  over  the  magnetic 
film,  performs  the  double  roles  of  providing  mode-conversion  compli¬ 
mentary  to  that  of  the  magnetic  film,  and  allowinp  the  possibility  of 
degenerate  mode  propagation  or  phase-matching.  Figure  5  illustrates 
that  the  bulk  wave  analogue  of  this  preferred  geometry  is  one  in 
which  the  anisotropic  and  optically  active  45°  rotators  are  sliced 
up  into  small  sections  which  are  then  interleaved. 

It  is  the  basic  design  of  Fig.  5  which  was  adopted  in  this 
program  as  the  best  candidate  for  experimental  study.  A  survey  of 
technically  proven  materials  reveals  that  the  best  choice  for 
experimentation  comprises  a  magnetic  garnet  film  grown  on  gadolinium 
gallium  garnet  with  LiI03  acting  as  the  anisotropic  super strate. 
Detailed  design  data  for  this  system,  based  on  the  theory  outlined 
in  the  next  section  is  presented  in  Section  4  of  this  report. 

3.  THEORETICAL  ANALYSIS  OF  OPTICAL  WAVEGUIDE  ISOLATOR 


In  this  section  we  outline  our  theoretical  analysis  of  an 
optical  waveguide  isolator/ gyrator  formed  by  a  magneto-optic  film 
sandwiched  between  an  isotropic  and  an  anisotropic  layer  as  shown  in 
Fig.  5.  The  mode- conversion  analysis  follows  the  ray-optical  approach 
of  Kogelnik  et  al^'  which  includes  the  Goos  Haenchen  effects.  That 
analysis  is  followed  by  a  brief  summary  of  how  the  anisotropic 
crystal  which  provides  mode-conversion  complimentary  to  the  magnetic- 
film  also  provides  for  degenerate  mode  propagation. 

3.1  Derivation  of  Coupled  Mode  Equations 

In  Fig.  6  we  depict  the  zig-zag  ray  path  of  a  TE  wave  in  the 
waveguide.  Mode  conversion  due  to  the  magnetic  film  will  occur  along 
the  paths  BC  and  DE  and  that  due  to  the  anisotropic  layer  will  occur 
upon  reflection  at  CD.  In  each  case  the  mode -conversion  may  be 
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described  in  terms  of  a  coupling  matrix  where  the  off-diagonal  terms 
represent  the  mode  conversion  coefficients.  In  general  we  write 


TE 

R, , 

R,  „ 

TE 

out 

11 

12 

in 

TM 

out 

R21 

22 

i 

i 

TM. 

in 

- 

L  J 

There  is  no  mode  conversion  at  the  substrate  boundary  because  the 
substrate  is  optically  isotropic.  The  mode  coupling  matrix  for 
reflection  at  the  substrate  boundary  therefore  has  zero  off-diagonal 
terms  and  is  given  by 


(3.2) 
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where  0  and  0  are  the  Goos-Haenchen  phase  shifts  associated  with 
total  internalreflection. 

The  anisotropic  top  layer  is  oriented  in  such  a  way  that  the 
permittivity  tensor  has  the  element  non- zero  when  expressed  in 
the  coordinate  system  shown  in  Fig  67  The  mode  coupling  matrix  for 
tlixs  reflection  is  given  by(2) 
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r,  8.  and  the  0's  have  been  related  to  material  parameters  by  Wang 
et  al^’  and,  for  convenience,  these  relations  are  tabulated  in 
Appendix  1. 


The  mode  coupling  matrix  applicable  to  a  single  passage  through 
the  magneto-optic  film  (either  B  to  C  or  D  to  E)  is  derived  in 
Appendix  2  and  is  given  by 
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where  &  =  - 
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In  adopting  the  analytical  approach  of  Kogelnik  et  al  we  use 
the  matrices  given  by  Eqs.  3.2,  3.3  and  3.4  to  write  down  the  TE  and 
TM  fields  at  each  end  of  a  complete  zig-zag  cycle.  These  are  then 
related  by  a  Taylor's  series  expansion  to  the  field  derivatives  and, 
after  some  algebra,  a  pair  of  coupled  differential  equations  emerge 
as 


f-B0 

AE=iKAM 

XD 
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Ag  and  ^  are  the  normalized  TE  and  TM  field  amplitudes,  (3  and 
are  the  propagation  constants  of  the  TE  and  TM  modes,  and  K  is  a 
coupling  coefficient  between  the  modes  given  by 
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bE  and  bM  are  the  round  trip  bounce  lengths  illustrated  in  Fig.  6. 

Equation  3.6  shows  us  how  non-reciprocal  mode  conversion  effects 
arise.  If  the  propagation  direction  is  reversed  with  respect  to  the 
magnetization  direction,  G,  and  hence  Y  will  change  sign.  If  we 
take  care  to  match  the  magnitude  of  the  two  terms  in  the  numerator  of 
Eq.  3.6  then  nor-reciprocal  mode  conversion  will  be  achieved;  K  being 
non-zero  in  the  forward  direction  (say)  but  zero  in  the  reverse 
direction.  The  interaction  length  for  complete  mode  conversion  may 
be  found  from  the  solutions  to  Eq.  3.5.  If  all  the  light  is  in  the 
TE  wave  at  Z  =  0  then^' 


vz> = 
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where  the  phase  mismatch  factor 
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The  maximum  converted  light  intensity  falls  from  100%  for  Q  =  0 
to  507o  for  Q  =  0.80.  Such  a  value  of  Q  allows  one  to  estimate  a 
practical  tolerance  on  the  phase-matching  between  TE  and  TM  modes. 
This  point  is  considered  further  in  Section  4. 


We  see  from  Eq.  3.7  that  complete  conversion  will  occur  in  a 
distance  -40pt  for  which  -L0ptK  =  it/2.  For  the  material  parameters 
given  in  Section  4  we  find  that  -i0pt  ~  2  mms. 

3.2  Phase-Matching  for  Degenerate  Modes. 


It  is  evident  from  Eq.  3.7  that  the  TE  and  TM  modes  between 
which  we  have  arranged  non-reciprocal  coupling  must  propagate 
degenerately  in  order  for  complete  TE  to  TM  (and  vice  versa) 
conversion  to  take  place.  The  diagonal  terms  of  the  mode- conversion 
matrices  given  by  Eq.  3.2,  3.3,  and  3.4  can  be  used  to  obtain  expres¬ 
sions  for  the  TE  and  TM  propagation  constants  j3g  and  f^.  In  order 
for  a  proper  waveguide  mode  to  exist,  the  phase  aelay  of  a  zig-zag 
wave  traversing  the  path  A-B-C-D-E  in  the  film  (Fig.  6  refers)  must 
be  the  same  as  that  of  an  evanescent  wave  travelling  in  the  substrate 
from  A  to  E.  After  some  algebra  this  requirement  reduces  mathematically 
to  the  following  equations. 
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where  n  is  the  effective  guide  index  P/k0  and  n  ,  n  ,  n  are  the 
refractive  indices  of  the  anisotropic  top-layer,  ^egisnerate  modes 
are  achieved  when  W_„  =  W^  for  a  special  value  of  n  .  Let  us 
examine  how  this  might  be  achieved  if  the  anisotropi§  crystal  is  a 
negative  uniaxial  crystal  with  its  optic  axis  in  the  film  plane, 
roughly  perpendicular  to  the  propagation  direction.  For  this 
configuration  n^  =  n^  =  nQ,  the  ordinary  refractive  index  and  n  =  n  , 


the  extraordinary  refractive  index.  Now  n^  >  n  for  any  waveguide 
and  this  gives  a  tendency  for  WT^  to  be  greater  thanW^p  for  any  value 
of  n  due  to  the  factor (n^/ng)<!f  in  the  second  arctan  term  in  Eq.  3.9. 
Sincl  in  a  negative  uniaxial  crystal  ne  <  n  we  may  find  a  value  of 
n  ,  close  to  nQ  for  which  (iip2/^2)2  (n  2  -°n02)  <  (n  2  -  ne2)  by  just 
tne  amount  to  counter  the  factor  (n^/n8y  yielding  the  result 
WfE  =  ^XM'  ^is  situation  is  illustrated  diagrammatically  in  Fig.  7. 

4.  DESIGN  DATA  FOR  LilO^  -  YIG  -  GGG  WAVEGUIDE  STRUCTURE 

In  order  to  attempt  an  experimental  demonstration  of  a  waveguide 
isolator  in  a  short  time  scale  (i.e.  this  reporting  period)  we  must 
select  our  materials  from  those  whose  growth  technology  has  been  well 
developed  for  other  reasons.  Magnetic  garnet  films,  grown  epitaxially 
on  to  gadolinium  gallium  garnet  (GGG),  have  received  a  lot  of  attention 
as  magnetic  bubble-domain  memory  media.  These  films  are  immediately 
suitable  as  optical  waveguides,  the  magnetic  films  having  a  higher 
refractive  index  than  GGG.  The  only  change  required  for  our  purpose 
is  to  move  the  magnetization  axis  from  normal  to  the  film-plane  to 
in  the  film-plane.  This  is  achieved  by  altering  the  film  composition 
slightly  to  change  the  lattice  parameter  mismatch  between  film  and 
substrate.  The  films  supplied  to  us  by  Dr.  R.  Henry  of  Rockwell 
International  (Anaheim)  have  the  composition  (LaY)^  Ga  33  Fe^  ^  0, 
and  are  grown  by  liquid  phase  epitaxy  on  to  <111>  oriented  G The 
transmission  spectrum  of  the  magnetic  garnets  demands  the  use  of  a 
1.06  pm  or,  for  even  smaller  absorption,  a  1.152  pm  laser  beam.  We 
use  a  cw  HeNe  laser  operating  at  1.152  pn  in  our  experiments.  Having 
settled  on  materials  for  the  film  and  substrate  the  anisotropic  top 
layer  has  to  be  chosen.  The  constraints  are  that  the  refractive 
indices  should  lie  within  a  range  (which  proves  to  be  small)  over 
which  mode -degeneracy  could  be  achieved.  This  phase-matching  require¬ 
ment  also  demands  that  the  material  shall  be  either  a  negative 
uniaxial  crystal,  or  a  biaxial  crystal  with  a  large  optic  angle.  It 
would  appear  that  lithium  iodate  LilOr  is  about  the  only  suitable 
material  with  a  developed  growth  technology. 

Final  design  data  may  be  obtained  from  the  theory  outlined  in 
section  3  once  the  refractive  indices  of  the  materials  are  known.  We 
know  from  section  3.2  that  the  optic  axis  of  LilO^  chould  be  set 
roughly  perpendicular  to  the  propagation  direction  in  order  to  most 
readily  achieve  degenerate  modes. 


We  first  assume  that  the  optic  axis  is  Indeed  at  90°  to  the 
propagation  direction  and  use  eqs.  3.8  and  3.9  to  calculate  the 
thickness  for  which  the  TE  and  TM  modes  are  degenerate.  In  bulk  YIG 
the  Faraday  effect  produces  about  280°/cm  optical  rotation  at  1.152  pm^^K 
This  corresponds  to  G  =  3.4  x  10"^  and  we  therefore  have  all  the 
information  needed  to  compute  the  magneto  optic  mode-coupling,  v,  given 
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by  eq.  3.4.  We  next  use  eq.  3.6  to  find  the  anisotropic  mode  coupling 
term,  S,  that  yields  zero  net  mode-conversion.  S  in  turn  is  related 
to  the  angle,  cc,  that  the  optic  axis  makes  with  the  propagation 
direction  via  the  equations  given  in  appendix  1.  It  turns  out  that 
a  ^  88  degrees  for  our  case  and  furthermore  the  phase-matching  is 
insensitive  to  a  as  shown  in  fig.  8.  This  fact  gives  us  the  experi¬ 
mentally  convenient  feature  that  a  can  be  readily  altered  by  suitably 
tilting  the  waveguide  assembly,  and  doing  this  will  "tune-in"  the 
correct  value  of  S  without  materially  altering  the  mode-degeneracy 
condition. 

The  required  design  data  is  summarized  in  Table  1.  In  this  table 
the  refractive  indices  of  the  garnet  film  and  substrate  were  measured 
in  the  laboratory  (see  sec.  5.11  and  those  for  the  LiI03  top-layer 
were  taken  from  the  li  terature*-5' .  All  other  data  emerges  from 
calculations  based  on  the  theory  given  in  section  3.  For  example  the 
tolerance  figure  on  the  design  film  thickness  of  1.395  im  is  computed 
from  the  allowable  phase-mismatch  introduced  by  letting  the  maximum 
converted  light  intensity  fall  from  1007,  to  50%.  It  should  be  noted 
that  a  maximum  of  507»  mode  conversion  would  still  allow  a  good 
demonstration  of  non-reciprocal  propagation  because,  in  the  reverse 
direction,  we  can  still  arrange  for  there  to  be  no  mode  coupling. 

5.  LABORATORY  TECHNIQUES 

Of  the  several  lab  techniques  and  practices  that  were  used  in 
the  course  of  this  study  we  select  the  following  three  items  as  being 
noteworthy.  The  first  describes  a  novel  prism  coupler  assembly 
which  allows  the  convenient  observation  of  optical  waveguiding  in 
small  samples  -  we  have  used  samples  of  less  than  3  mns  long  in  any 
dimension.  The  second  item  is  concerned  with  how  a  knowledge  of 
coupling  angles  to  excite  waveguide  modes  can  lead  to  accurate  film 
thickness  and  refractive  index  determinations.  Lastly  we  report  the 
results  of  an  experiment  in  which  the  thickness  of  a  garnet  waveguide 
film  was  reduced,  in  a  controlled  way,  to  a  specified  value  using  a 
sputter-etch  technique. 

5.1.  Prism  Coupling  to  Small  Pieces  of  Optical  Waveguides 

Rutile  has  a  sufficiently  large  refractive  index  to  couple  to 
all  propagating  modes  in  any  garnet  film  waveguide.  Rutile  prisms 
available  to  us  were  of  the  45°/45°/90°  type  with  a  1cm  square 
aperture.  Such  prisms  would  normally  be  larger  than  the  waveguide 
samples  and  we  therefore  designed  a  mounting  device  in  which  an  input 
and  an  output  prism  formed  a  planar  mounting  platform  for  the  wave¬ 
guide  films.  Small  plastic  rods  transmitted  the  necessary  coupling 
pressure  to  the  right-angle  corners  of  the  prisms.  A  photograph  of 
the  assembly  and  an  exploded  view  is  shown  in  fig.  9,  and  a  sketch 
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of  the  metallic  components  is  given  in  fig.  10.  The  planar  mounting 
surface  was  obtained  by  glueing  the  prisms  to  their  brass  holders  and 
then  rubbing  the  latter  on  emery  cloth  until  observation  under  a  high 
power  microscope  indicated  Lhat  both  prism  top- surfaces  were  in  the 
same  plane.  The  whole  assembly  depicted  in  fig.  9  was  mounted  on  a 
goniometer  structure  whose  rjtation  axis  could  >-e  adjusted  to  coincide 
with  the  right  angle  corner  of  the  input  prism.  In  this  way  large 
changes  in  incidence  angle  could  be  made  without  the  need  to  reposition 
the  beam  to  the  prism  corner. 

When  coupling  to  waveguides  that  are  formed  on  thin  substrates 
confusion  can  arise  in  identifying  coupled-out  waveguide  modes  from 
light  that  gets  into  the  substrate  at  the  first  prism  and  out  again 
at  the  output  prism.  The  substrate  radiation  is  exited  by  the 
incident  beam  being  scattered  into  a  large  angular  spread  by  the 
finite  size  of  the  right-angle  corner  of  the  input  prism.  We  mount 
a  screen  on  top  of  the  goniometer  which  therefore  rotates  with  the 
prism/waveguide  assembly  as  the  incidence  angle  is  varied.  Waveguide 
modes  impinge  on  this  screen  at  fixed  locations  whereas  spots  of 
radiation  from  the  substrate  move  as  the  incidence  angle  is  adjusted. 
Spots  of  light  originating  from  waveguide  modes  are  therefore  readily 
distinguishable.  Since  we  use  a  1.152  jan  laser  beam  all  our  observa¬ 
tions  were  made  with  the  aid  of  an  image  converter  tube  with  an  SI 
photocathode. 

5.2  Optical  Waveguide  Technique  to  Determine  Film  Thickness  and  Index. 

In  order  that  we  might  accurately  specify  the  waveguide  film 
thickness  we  must  have  good  data  on  refractive  indices,  etc. 

Fortunately  the  waveguides  themselves  offer  a  good  technique  for 
refractive  index  and  film  thickness  determination.  The  procedure  is 
to  accurately  measure  angles  at  which  waveguide  modes  are  exited. 
Knowing  the  refractive  indices  and  prism-angle  of  the  coupling  prism 
we  can  calculate  the  propagation  constants  of  these  modes  and  then  fit 
this  data  to  the  mode  dispersion  equations  (eqs.  3.8,  3.9)  using  a 
least  squares  algorithm.  Our  prism  data,  obtained  using  a  spectrometer 
operating  at  .6328  and  1.152  |um  is  summarized  in  table  2.  The  angles 
of  incidence,  9^  at  which  waveguide  modes  are  exited  are  related  to  the 
mode  propagation  constant  g  by  the  equation 


g/kQ  =  n^  sin(a  +  arcsin(sin9i/np) ) 


(5.1) 


where  n  is  the  prism  refractive  index  and  a.  is  the 

D 

incident  and  waveguide  faces  of  the  prism. 


angle  between  the 
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We  programmed  a  Hewlett  Packard  Model  20  calculator  to  perform 
the  function  fitting.  The  program  operates  on  the  principle  of 
selecting  the  value  of  film  refractive  index  that  gives  the  smallest 
standard  deviation  of  film  thickness  determinations  from  each  measured 
value  of  f.  Usually  the  substrate  refractive  index  is  entered  as  a 
given  value  although  it  is  possible  to  also  vary  that  parameter  to 
find  the  best  fit.  Obviously  we  have  to  record  coupling  angles  for 
at  least  as  many  parameters  that  we  wish  values  for  (i.*.  .  n^  and 
W,  or  nf,  W  and  ng). 

During  our  work  we  found  that  mea0uring  several  modes  (3  to  6) 
to  about  ±  5  arc  minutes  accuracy  would  yield  data  with  a  standard 
deviation  of  about  ±  .003  gm  on  thickness  and  ±  .0004  on  refractive 
indices.  Data  for  typical  experiment  is  given  in  table  3.  The 
goniometer  reading  at  normal  incidence,  taken  when  the  incident  lase1' 
beam  was  reflected  back  along  its  path,  served  to  calibrate  the 
goniometer.  Given  the  apparatus  we  find  that  a  new  film  can  be 
accurately  characterized  in  less  than  half  an  hour,  most  time  being 
spent  in  accurately  fixing  the  coupling  angles.  We  feel  therefore 
that  the  technique  described  above  offers  a  fast  and  accurate  method 
to  determine  refractive  indices  and  thicknesses  of  a  rather  special 
type  of  dielectric  thin  film,  namely  a  film  of  up  to  10  ym  thick 
grown  or  deposited  on  to  a  lower  refractive  index  substrate. 

5 . 3  Sputter  Etching  Experiment  to  Reduce  Film  Thickness. 

The  final  design  data  for  the  waveguide  isolator  experiment 
calls  for  a  film  thickness  of  1.395  ±  .02  4m.  This  figure  depends 
upon  the  refractive  index  of  the  film.  One  of  our  samples  (# 6/2A/8) 
came  to  us  from  Rockwell  International  at  1.85  4m.  We  decided  to 
remov-s  material  from  this  film  by  sputter  etching.  Our  equipment  was 
a  VEEC0  Model  716  evaporating  plant  modified  for  RF  sputtering  at 
13.56  MHz.  The  sputtering  atmosphere  was  Argon  gas  at  about  33  4m 
pressure.  The  anode  to  cathode  separation  was  15  mms.  We  chose  to 
run  at  50  watts  RF  power  for  a  standard  15  minute  period.  After 
each  period  the  sample  was  annealed  in  air  at  800° C  for  30  minutes. 

(We  found  that  the  sample  was  reduced  somewhat  in  the  sputtering 
apparatus  due  to  heating  ir  1  low  pressure  environment.) 

We  repeated  the  above  procedure  5  times  to  reduce  the  sample 
thickness  from  1.847  4m  to  1.382  4m.  Th'  measured  data  is  summarized 
in  table  4.  Notice  that  the  refractive  icdex  remained  unchanged 
throughout  the  program.  Figure  11  shows  the  sample  thickness  as  a 
function  of  etch  time.  We  have  no  explanation  for  the  change  in 
etch  rate  halfway  through  the  experiment.  Every  effort  was  made  to 
keep  each  sputter  etch  run  the  same  in  terms  of  gas  pressure,  RF 
power,  etc.  Surface  roughness  did  not  deteriorate  markedly  during  the 
etching  program. 
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6 .  EXPERIMENTS  ON  "ONE-WAY"  MODE  CONVERSION 


Our  experimental  program  to  study  non-reciprocal  effects  in 
optical  waveguides  is  centered  or  a  GGG-YlG-LilO^  waveguide  of 
appropriate  dimensions  and  a  pair  of  rutile  right-angle  prisms 
acting  as  input/ output  couplers  as  depicted  in  fig.  12.  We  apply  a 
magnetic  bias  field  parallel  to  the  propagation  direction,  couple 
into  the  TE  mode  and  look  at  the  amount  of  TM  light  that  emerges. 
Non-reciprocal  effects  will  be  evident  if  the  amount  of  mode-conversion 
changes  when  the  bias  field  is  reduced.  The  LilO^  orientation  could 
be  altered  to  minimize  the  smaller  amount  of  conversion. 

We  mounted  our  waveguide  on  the  rig  described  in  section  5.1. 

The  whole  assembly  was  made  of  non-ferrous  materials  and  so  we  were 
able  to  apply  bias  fields  of  about  100  Gauss  with  a  suitably  located 
electromagnet.  A  photograph  of  the  assembled  apparatus  is  given  in 
fig.  13.  An  index-matching  liquid  under  the  rutile  prisms  minimizes 
mismatch  reflections  between  the  prism/YlG/GGG  and  LHO3 -YlG-GGG 
regions. 

Unfortunately  our  experiments  h*ve  not  been  successful  in  demon¬ 
strating  non-reciprocal  effects.  We  attribute  this  to  our  failure  to 
get  I ilO^  and  garnet  film  in  close  enough  contact.  We  have  seen  small 
amounts  of  mode  conversion^),  insensitive  to  bias  field  reversal,  but 
the  TE  and  TM  modes  are  not  degenerate,  indicating  that  the  anisotropy 
of  the  LiI03  is  not  playing  its  part. 

We  have  studied  the  mode  dispersion  waves  for  a  4  layer  waveguide 
where  va  have  introduced  a  thin  gap  of  thickness  t  between  the  LiI03 
and  the  garnet  film.  In  fig.  14  we  plot  TE  and  TM  curves  for  several 
values  of  t.  We  see  that  the  gap  is  a  critical  factor  in  the  location 
of  the  degenerate  mode  point.  What  happens  is  that  the  waveguide  mode 
field  decays  into  the  gap  and  consequently  the  LiI03  is  quickly 
screened  from  the  mode  It  only  takes  a  t  =  0.1  |»  to  remove  any 
prospects  of  mode  degeneracy.  A  gap  of  less  than  0.1  fan  infers  a 
good  optical  contact.  Wo.  could  not  optically  contact  our  samples;  a 
better  optical  finish  will  have  to  be  obtained. 

Alternatively  we  could  place  a  high-index  liquid  layer  between 
the  LiI03  and  the  garnet  film.  If  the  r«._ractive  index  of  this  liquid 
were  greater  than  the  effective  guide  index  of  the  waveguide  modes  of 
interest  then  we  would  get  a  sinusoidal  field  pattern  in  the  liquid 
layer  instead  of  an  exponential  decay  and  the  LiI03  would  continue  to 
be  effective.  We  have  such  a  liquid  to  hand;  it  is  a  Cargille  Labora¬ 
tories  index  matching  liquid  from  their  EH  series.  Its  refractive 
index  for  the  sodium  D  line  is  2.04,  we  measure  n  =  1.973  and  a  =  0.9  cm 
at  1.152  jim.  Unfortunately  the  material  contains  arsenic  tribromide 
which  reacts  with  the  LiI03  and  so  we  cannot  use  it.  We  have  explored 
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the  location  of  the  degenerate  mode  operation  points  as  a  function  of 
the  thickness  of  a  liquid  layer  of  index  1.973.  We  see  that  the 
possibilities  for  degenerate  modes  persist  to  quite  thick  films  as 
shown  in  fig.  15.  Furthermore  the  phase-matching  condition  is 
insensitive  to  the  exact  liquid  layer  thickness  if  it  is  the  range 
0.5  to  1.0  jum. 


SUMMARY 

During  this  reporting  period  (1  December  1972  to  31  January  1974) 
we  have  identiried  the  requirements  for  a  magneto- optic  isolator  in  an 
optical  dielectric  waveguide.  These  are  that  the  provision  of  a  magneto¬ 
optic  medium  is  necessary  but  not  sufficient.  We  must  also  provide 
an  anisotropic  medium  to  give  us  the  capability  of  achieving  degenerate 
modes  in  the  waveguide,  and  to  provide  a  mode  coupling  coefficient 
complimentary  to  the  magnetic  one.  A  periodic  structure  like  Tien's 
serpentine  circuit^8'  cannot  be  used  to  make  up  for  a  mode  mismatch 
since  any  such  structure  destroys  the  unique  symmetry  of  the  longi¬ 
tudinal  magneto-optic  effect  which  produces  non-reciprocal  phenomena. 

We  have  carried  out  a  detailed  theoretical  study  of  the  problem  and 
are  confident  that  our  design  is  the  best  choice  based  on  currently 
available  materials .  We  use  a  layered  structure  in  which  an  iron 
garnet  film  is  grown  epitaxially  on  to  GGG  and  a  LilO-  cyrstal  is 
optically  contacted  to  the  film  to  form  a  superstrateT  We  have  been 
successful  in  getting  all  the  materials  to  the  requived  dimensions  etc. 
and  have  demonstrated  our  ability  to  measure  critical  parameters  like 
film  thickness  to  adequate  accuracy. 

So  far  we  have  failed  to  demonstrate  non- reciprocal  effects. 

We  are  confident  that  this  is  due  to  our  inability  to  optically  contact 
the  LiI0_  to  the  garnet  film.  The  available  samples  had  somewhat 
relaxed  tolerances  on  optical  flatness  etc.  New  samples,  with  better 
optical  finish  will  become  available  and  we  look  forward  to  better 
success  in  optically  contacting  these  satisfactorily. 

An  alternative  avenue  for  further  research  would  be  to  explore 
the  possibilities  of  non-isotropic  magneto-optic  materials  in  which 
the  magnetic  and  anisotropic  mode  coupling  would  take  place  in  the 
same  film.  At  this  writing  we  have  not  begun  to  analyze  such  systems 
in  any  detail.  A  possible  material  class  would  be  the  orthoferrites. 
Study  would  be  required  to  determine  whether  the  phase-matching 
possibiliti  ^  could  lead  to  experimentally  attractive  situations  like 
that  found  L.  the  system  described  in  this  report.  We  refer  here  to 
the  happy  situation  illustrated  in  fig.  8  where  orientation  adjustments 
change  the  mode  coupling  coefficient  by  large  amounts  whereas  the 
phase-matching  is  hardly  affected.  Preliminary  investigations  suggest 
that  this  type  of  behavior  does  not  normally  occur  in  waveguides  with 
anisotropic  films. 
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APPENDIX  1 


MODE  CONVERSION  MATRI/  FOR  REFLECTION  AT  ANISOTROPIC  DIELECTRIC 

BOUNDARY 


Wang  et  alvw  have  derived  expressions  relating  the  reflection  and 
mode  conversion  that  takes  place  between  TE  and  TM  waves  upon  total 
internal  reflection  at  a  boundary  between  an  isotropic  medium  of  high 
refractive  index  and  an  optically  uniaxial  medium  whose  optic  axis 
lies  in  the  boundary  plane  For  convenience  their  expressions  are 
copieo  out  below. 
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In  the  above  expression  and  p2  are  the  decay  constants  of  the 
evanescent  waves  into  the  anisotropic  layer,  given  by 
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n  is  the  effective  guide  index  £/K  and  the  K^j  are  the  permittivity 
tlnsor  coefficients  for  the  uniaxia?  layer.  If  the  optic  axis  makes 
an  angle  of  a  with  the  propagation  direction  we  find  that 
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r  ,  r  and  r  are,  in  general,  complex  quantities.  In  section 
3  of  tflfs  rlport  they  are  referred  to  in  a  modules  and  phase  angle  form. 
Wang  et  al^2l  also  show  that  we  may  write 
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For  a  given  configuration  r,  S,  0e  and  0^  are  readily  obtained  by 
numerical  computation  from  eqs.  1.2. 
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APPENDIX  2 

MODE  CONVERSION  MATRIX  FOR  CUBIC  MAGNETO-OPTIC  WAVEGUIDING  FILM 

The  amount  of  TE  to  TM  mode  conversion  (and  vice  versa)  that 
takes  place  during  a  single  traverse  of  the  "zig-zag"  plane  wave  that 
makes  up  the  waveguide  mode  may  be  described  in  terms  of  a  matrix  R 
where 
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In  the  following  analysis  we  work  out  expressions  for  the  coefficients 
of  R  in  terms  of  the  properties  of  the  magneto-optic  film,  which  we 
assume  to  be  optically  isotropic. 

In  general  the  dielectric  properties  of  a  magneto- optic  medium 
may  be  described  by  a  gyration  vector  which  features  in  the  equation 
relating  D  and  E  as 


D=€  E  +  i  E  x  £ 


(2.2) 


For  an  isotropic  medium  the  dependence  of  £  on  the  external  field  is 
one  of  simple  proportionality,  J  too  reduces  to  a  scalar  n2,  where  n 
is  the  refractive  index.  We  deTine  a  right-hand  set  of  coordinate 
axes  x^(i=l,2,3)  with  X3  parallel  to  the  waveguide  propagation  direction, 
X2  also  in  the  film  plane  and  xi  directed  normal  to  the  film  form 
substrate  to  top  layer  (fig.  6  refers).  The  dielectric  constant 
tensor  referred  to  these  axes  is  therefore  given  by 


3 

=  Z 


j-1 


Ku  ej 


+  i 


3  3 

Z  Z 

j=l  k=l 


6ijkEj8K 


(2.3) 


where  e^jjc  is  an  antisymmetrical  unit  tensor. 

If  we  assume  that  the  magnetization  is  directed  along  the 
propagation  direction,  x^,  then  g.  =  =  0  and  g^  =  G(say).  The 

dielectric  tensor  then  takes  on  the  special  form 
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(2.4) 


It  will  be  convenient  to  transform  the  problem  to  new  axes  (y^) 
where  >3  is  directed  along  the  "zig-zag"  ray  path  in  the  film.  If 
y„  is  set  parallel  to  then  the  TE  and  TM  electric  fields  will  be 
directed  along  y  and  y^  respectively.  The  transformation  matrix 
relating  (x  )  and  ( y )  if  given  by  y£  =  ajaxj.  The  ±  superscript 
refers  to  tne  "up"  and  "down"  rays  respectively.  If  8  is  the  angle 
of  incidence  of  the  "zig-zag"  wave  (fig.  6)  then  we  write 


co:  (n/2  -  0) 


±sin(jt/2-  9) 


0  T  sin(it/2  -  9) 

1  0 

0  cas(r t/2  -  8) 


(2.5) 


We  can  now  refer  tha  dielectric  tensor  of  the  magnetic  film  to  the 
axes  (y^)  and  find  that 


n2  iGsin9 

-Gsin9  n2 


0 

TiG  cos  8  . 


(2.6) 


0 


±iG  cacQ 


2 

n 


If  we  assume  an  e*-^  "  —  ^  dependence  of  the  fields  the  electro¬ 
magnetic  wave  equations  reduces  to 


k2  E  -  k(k  •  E)  =  kQ2  Ky*  E  .  (2.7) 

Let  us  assume  a  solution  of  the  form  e^^1^  where,  for  case  of 

indexing  we  replace  y-  by  Z.,y.  by  y?and  y1  by  x.  We  therefore  obtain 
from  eq.  2.7  that 
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V2(i  Ex  +  i 


n2  iGsinO 

fi  E  ' 

2 

-iGsin0  n 

—  X 

i  e 

-i  y 

L  J 

and  hence  y2  =  (n2  ±  G  sin0)  and  E  =  ±  iE  . 

x  y 


(2.8) 


(2.9) 


Equation  2.9  shows  us  the  well  known  result  that  the  eigenmodes  of 
propagation  are  circularly  polarized  with  propagation  constants 


yr  =  (n2  ±  Gsin8)^. 


(2.10) 


Notice  that  the  same  result  is  obtained  for  both  the  "up"  and 
"down"  rays.  In  order  to  work  out  the  mode -conversion  that  takes  place 
in  one  transit  (up  or  down)  of  the  zig-zag  path  we  decompose  each  of 
the  TE  and  TM  polarizations  into  two  counterrotating  circularly 
polarized  modes.  We  let  these  propagate  along  the  ray  path  length 
(U leas'd  and  then  reform  the  TE  and  TM  polarizations.  Thus 


E  (0)  -iknyj  F,  (0)  -ik  Y  1  E  (0)  -ikQV  r  E  (0)  -ikoY 
F  (z)  =  _i_  e  +  -  e  +  i  -  e  -  i-* -  e 


(2.11) 


E  (0)  -ik  yj  E  (0)  -ik  E  (0)  -ik  Y*  E  (0)  -ik  y 


Ey(z)  =  -i  -f 


l 
o  z 


O'Z.  .  X  o'z+^__e  O  z+  - e  ~  (2tl2) 


+  i 


At  z=0  we  find  that  E  (z=0)  -  E  (0)  and  E  (z-0)  =  E  (0).  Equations  10  and 
11  can  be  factorized  by  noting,  for  example,  that 


r  -i(yr-ySkoz  i(yr+y^)kQz 
■ikoy  z  =  - 2  -  2 


Thus 


•ik0(Yr+Yi)z 


r»„<Yr-vVl  [k0(Yr 

L - 2 - r  Ey(0'  SlnL - 2 


Ex<2)  *  e 


E  (0)  COS 
x 


-ySz 


(2.13) 
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5y(z)  =  e 


■ikQ(v  +Y  )z 


■E  (0)  sin 
x 


[k0<Yr-/)zl  rko(Vr-*. 

[—2 - J  +  V°>  <*[-*— 2 - J 


(2.14) 


From  eq.  1 . 10  we  find 


r L 
V  +Y 


n  and 


Y  -Y  _  2GsiT>9 


and,  since  we 


1  X  UU1  «  XV  VYC  X  XUU  -  ***»'*  2  f  OXUVC  WC 

can  identify  TE  modes  with  E  and  TM  modes  witR  E  ,  we  find  that  the 
matrix  defined  in  eq.  2.1  can  be  written  from  eq.  2.13  and  2.14  as 


ri  4' 


DOWN]  i& 

J"  6 


cos  Y  - 

sin  Y 


sin  Y* 

OJY 


2.15 


where  &  = 
matrix. 


-k  r  V 
_ o _ 

casQ 


and  Y 


-kQ  g  w  tan  9 
n 


This  is  the  desired 
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TABLE  1 


DESIGN  DATA  FOR  LilO-- YIG-GGG  WAVEGUIDE 
ISOLATOR  FOR  1.152  \i m  OPERATION 

Substi ate 

Single  crystal  Gd^  Ga^  °12 

Film  grown  epitaxially  on  to  <111>  surface 

n  =  1.945  ±  .001 
s 

Magneto- optic  Film 

(LaY)3  Ga  33  Fe^O^  grown  by  LPE 
nf  =  2.1783  ±  .0004 
W  =  1.395  ±  .02  tun 

G  =  3.4  x  10-4  (=  280° /cm  Faraday  rotation) 

Y  =  2.43  x  10  3  rad 

m  =  2  modes  are  degenerate,  8/kQ  =  1.9568,  6  =  63.9° 
Top  Layer 

Single  crystal  LUO3  (optically  contacted) 

n  -  1.858,  n  =  1.716 
u  o  e 

Optic  axis  about  88°  from  2  in  yz  plane 

Length  for  complete  mode  conversion  =“  2  mms 

0  -  0  =  89.0°  -  82.8°  =  6.2° 

ee  vnm 

S  =  0.0049 

Waveguide  decay  constants  into  top  layer 


(a)  ■.or  TE  mode  1.08  pn 

(b)  for  TM  mode  1.70  tan . 
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TABLE  2 


REFRACTIVE  INDICES  OF  SYNTHETIC  RUTILE  AT  1.152  jim  and  .6328  jot 


WAVELENGTH 

n 

n 

“O 

— e 

.6328  jan 

2.5833  ±  .0003 

2.8647  • 

.  0003 

1. 152  jim 

2.4703  *  .0003 

2.7254  • 

.0001 

Note:  Measurements  were  made  at  24°C. 


TABLE  3 

Typical  data  taken  to  determine  film  thickness  and  refractive 
index  at  1.152  jt.  /  . 

Prism  angle  -  45°  29'.  Substrate  refractive  index  -  1.945. 

Sample  No.  6/2A/7.  Composition  (YLa).jGa  33Fe4  67°12' 

Scale  reading  at  normal  incidence  75°  48' . 


Scale  Reading 

56°40'±1' 

62°09'±2' 

70°31'±5' 

34°57'±1' 

43°39'±2' 

55° 16' ±5' 

Least  squares  fitting  yeild 


Mode  No. 

9 

B/K 

TE 

19.13° 

2.159 

TE? 

13.65° 

2.101 

TE^ 

5.28° 

2.007 

TM 

40.85° 

2.157 

TM? 

32.15° 

2.093 

tm.; 

20.53° 

1.990 

thickness  = 

1.715  ± 

.008  jim 

index  = 

2.1783  + 

.0006 
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TABLE  4 


Data  from  sputter  etching  experiment. 
(YLa)3Ga_33Pe4_67°12. 


Sample  No.  6/2A/8,  Composition 


UNIT 

ETCH  TIME  (MINS) 
TOTAL 

PRISM 

TEq 

COUPLING 

tel 

ANCLES  6 

te2 

W 

(ton) 

n 

0 

0 

19.40 

14.55 

6.95 

1.847 

2.178 

15 

15 

19.2/ 

14.07 

6.03 

1.771 

2.179 

15 

30 

19 . 0? 

13.47 

4.87 

1.688 

2.178 

15 

45 

18.09 

13.00 

3.80 

1.625 

2.178 

15 

60 

18.63 

11.95 

1.87 

1.502 

2.178 

15 

75 

18.43 

10.80 

0.15 

1.382 

2.179 
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FARADAY  EFFECT 


Figure 


Figure 


METAL  FILM 


Figure 


TOP  LAYER  (ANISOTROPIC) 
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Figure  14 


